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The epoch when the Universe had a temperature higher than a GeV is long before any time 
at which we have reliable observations constraining the cosmological evolution. For example, the 
occurrence of a second burst of inflation (sometimes called thermal inflation) at a lower energy 
scale than standard inflation, or a short epoch of early matter domination, cannot be ruled out by 
present cosmological data. The cosmological stochastic gravitational wave background, on scales 
accessible to interferometer detection, is sensitive to non-standard cosmologies of this type. We 
^\ ' consider the implications of such alternative models both for ground-based experiments such as 

0^ , LIGO and space-based proposals such as LISA. We show that a second burst of inflation leads to a 

■ scale-dependent reduction in the spectrum. Applied to conventional inflation, this further reduces an 

already disappointingly low signal. In the pre big bang scenario, where a much more potent signal 
is possible, the amplitude is reduced but the background remains observable by LISA in certain 
parameter space regions. In each case, a second epoch of inflation induces oscillatory features into 
the spectrum in a manner analogous to the acoustic peaks in the density perturbation spectrum. On 
LIGO scales, perturbations can only survive through thermal inflation with detectable amplitudes if 
CO ' their amplitudes were at one time so large that linear perturbation theory is inadequate. Although 

for an epoch of early matter domination the reduction in the expected signal is not as large as the 
one caused by a second burst of inflation, the detection in the context of the pre big bang scenario 
may not be possible since the spectrum peaks around the LIGO frequency window and for lower 



frequencies behaves as f 3 . 



PACS numbers: 04.30.-w, 98.80.Cq Preprint |gr"-qc/981104c| 



CN 
> 
O 

o 

00 
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A stochastic gravitational wave background is potentially a powerful source of cosmological information |1[ , since 
gravitational waves decoupled from the matter in the Universe at an extremely early stage. Various ways exist in 
(sj). which such a background might be generated, including cosmological inflation, collision of bubbles during a phase 
' transition, and motions of topological defects. We will be concentrating on the first of these. 

There are currently four large-scale ground-based interferometers under construction, namely LIGO, VIRGO, GEO 
and TAMA, and a proposal for a space-based detector LISA is also under consideration; see Ref. Jjj for reviews of these 
experiments. Unlike measures of density perturbations, these experiments probe frequencies which are high enough 
that the gravitational waves have had a wavelength less than the Hubble radius since very early in the history of the 
Universe. Assuming the standard cosmology, where the Universe was radiation dominated from an early stage, giving 
way to matter domination relatively recently, the time when a comoving frequency /* equalled the Hubble scale is 
given by 

/* i?*0* ^ z i/2 

fo Hoao T cq eq 

where f = a H n = 3ft- x 10~ 18 Hz is the mode that is just re-entering the Hubble radius today, T is the temperature 
and z oq = 24 000 fl^h 2 the redshift of matter-radiation equality. Here flo and h are the density parameter and 
Hubble constant in the usual units. Comoving units have been normalized to match present physical units. Since 
T cq = 24 000 tt h 2 T « leV, we have 



eq 



100 Hz 10 9 GeV 



(2) 



where 100 Hz is the typical frequency band probed by ground-based interferometers. 

The epoch when the Universe had a temperature higher than a GeV is long before any time at which we have reliable 
observations constraining the cosmological evolution. Because interesting wavelengths have been less than the Hubble 
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radius since then, the amplitude of the stochastic background will be sensitive to whatever the actual cosmological 
evolution is. A different cosmological model will have two principal effects: first, it will modify the relationship 
between scales given by Eq. (0) for the standard cosmology, and secondly it can alter the amplitude. Only those 
scales which are larger than the Hubble radius will have their amplitude unchanged in alternative cosmologies. 

To highlight this effect, we shall consider two alternative scenarios to the standard cosmology. In the first we assume 
a second short burst of inflation at a lower energy scale than standard inflation, often called thermal inflation [||. 
Such a short stage of inflation may be desirable to get rid of unwanted relic particles which could otherwise have 
dramatic consequences on the subsequent evolution of the Universe |2j . This second burst of inflation has in general no 
connection to the usual inflationary epoch which is assumed responsible for the generation of the density perturbations 
and gravitational waves leading to structure formation and microwave anisotropies. It is normally imagined to begin 
at a temperature around 10 7 GeV, and continue until about 10 3 GeV (the supersymmetry scale), giving a total of 
In 10 4 ~ 10 e-foldings of inflation. If this second burst of inflation occurs after the scales of interest (i.e. the scales 
accessible to LIGO and LISA) are within the horizon, then the stochastic background is altered. 

The second alternative is the recently discussed possibility that the Universe becomes temporarily dominated by 
some long-lived massive particle (for example the moduli fields of string theory), inducing an early period of matter 
domination which ends when the particle finally decays to restore the usual radiation-dominated Universe (see Ref. || 
and references therein). Another possible origin for a short period of early matter domination lies with the QCD 
phase transition, as was recently discussed by Schwarz jl] . 

To make the discussion concrete, we require a model for the initial gravitational wave spectrum. We consider two 
somewhat different models for the actual origin of the stochastic background. The first is conventional inflation, 
which is the standard case but which gives a disappointingly low signal even in the standard cosmology. The second 
is an unusual variant on inflation known as the pre big bang (PBB) scenario |5|J1. This is a much more speculative 
idea, but potentially of greater interest as it predicts a spectrum which rises sharply towards small scales , the 
opposite of the usual inflationary behaviour, thus creating a realistic prospect of direct detection of the gravitational 
wave background ||. 



II. COSMOLOGY WITH THERMAL INFLATION 



The first modification to the standard Big Bang model we consider is a second epoch of inflation occurring after 
the beginning of the radiation epoch. We assume that the original epoch of inflation (either conventional or pre big 
bang) responsible for generating the gravitational wave background ends at conformal time r r i (from now on we will 
always use conformal time), at which point the radiation era begins. At time r i2 a second era of inflation starts which 
ends at r r2 ; we will model this era as exponential inflation. After r r2 the Universe again becomes radiation dominated 
and the evolution of the Universe continues according the standard model with an era of matter domination starting 
at r m which lasts up to the present. The scale factor after the original period of inflation is given by 



' a 2 (r 2 + t) 
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r r i < r < T i2 



03 (73 - r) 1 ; r i2 < r < T r2 
a 4 r ; r r2 < r < r cq 
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where the constants are determined by requiring the continuity of a and da/ dr at the transition between two epochs. 
For the moment we only want to make an approximate estimate of the modifications introduced in the spectrum of 
the stochastic gravitational wave background and we will not worry about the exact value of these constants. We 
parametrize the thermal inflation by the energy scale at which it starts, pi 2 , and by the total expansion it gives rise 
to, i?i = a r2 /a i2 . 

In this model we assume that the first epoch of radiation domination happens at a temperature of order I0 16 GeV, 
while the short period of thermal inflation begins around 10 8 GeV. We also assume that during thermal inflation the 
Universe expands by a factor of I0 4 (approximately 10 e-foldings of inflation). This is the minimal amount of inflation 
which solves the moduli problem J2|. 

Normalizing the scale factor at the present as ao = 1, we can write all the relevant parameters in our model in 
terms of H u Hj 2 (the Hubble parameters at the end of the first period of inflation and at beginning of the second 
stage of inflation), i?i, the present Hubble parameter Hq and the redshift at the matter-radiation equality z cq . In 
terms of these parameters we can write the Hubble parameter and time at the matter-radiation equality 
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H cq = (1 + z cq )^ 2 H , r cq = (1 + z cq y 1/2 H^ 1 , (4) 
the Hubble parameter, redshift and time at the end of the second epoch of inflation 



,1/4 (Hi2\ 1/2 _ (1 + Z cq ) 1/4 



H l2 =H i2 , 1 + z r2 = (1 + z cq ) 1/4 , r ti = yjT( u , , (5) 



as well as the redshift and time at the beginning of the second period of inflation 

(Ha\ 1/2 Ri (1 + 

\H ) ' Ti2 ~ (H^Ho) 1 / 2 



l + z i2 = (1 + Zcq) i?i ( — ) , T i2 « 77TT7TT72-' ( 6 ) 



and at the end of the first epoch of inflation 
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We now describe the physical effects through which thermal inflation modifies the spectrum. The technical details 
of the calculation are quite complex, and we describe them fully in the Appendix; here in the main text we will restrict 
ourselves to describing the relevant physics and exhibiting the results. 

On short and long scales things are relatively simple. Modes always above the horizon cannot change their amplitude, 
but there is an apparent effect caused by the change in the correspondence between comoving scales (i.e. physical scales 
during inflation) and present physical scales; the extra expansion during thermal inflation shifts the spectrum to long 
wavelengths. For a flat spectrum this has no effect, but if the spectrum is tilted this alters the amplitude on a given 
present physical scale. This effect may have observational consequences in the contribution of tensor perturbations 
to the anisotropy of the cosmic microwave background radiation as we will discuss in the following subsection. 

A related effect is a redshift in the high-frequency cutoff of the spectrum corresponding to the horizon size at the 



end of the main period of inflation. From Eq. (A21) we obtain 



_ (H rl H Q f 2 

/rl - (1 + ^)^ ■ (8) 

Very short scale modes, which are always inside the horizon when the alternate behaviour is occurring, simply 
suffer an additional redshift. Modes inside the horizon behave as radiation, with the energy density decreasing as 
a~ 4 . If thermal inflation occurs it means there is more total expansion between the initial generation of the modes 
and the present, and so their amplitude is suppressed. For example, with perfect exponential thermal inflation, this 
suppression will be Rr 4 ~ I0~ 16 . Had we assumed an epoch of power-law thermal inflation with aocr 9 (q > 1), 

/ r i would be redshifted by j^ 1+q ^ 2q , Recalling that during an epoch of thermal inflation the total energy density 

redshifts as a 2 ^- 1 ^ p ^ 1 \ the suppression in fi gw in this case would be by a factor Pgw/p = -R; 2 ^ 1+1 ^ q \ 

The most interesting regime is the intermediate one. For waves with frequencies in the interval f l2 < / < f r 2, a 
curious phenomenon happens: these waves enter the Hubble radius during the first epoch of radiation but are again 
pushed outside the horizon during the subsequent epoch of inflation. While they are inside the horizon their behaviour 
is oscillatory. Hence, when they leave the horizon, their phases will have a strong dependence on frequency and this 
will leave an oscillatory imprint in the spectrum. This is strongly reminiscent of the acoustic peak structure in the 
density perturbation spectrum which leads to oscillations in the radiation power spectrum, and is arising similarly 
from the dominance of growing mode perturbations at horizon entry. Similar oscillations have also been found in the 
primordial density perturbation spectrum in the case of double inflation . 

To make our discussion concrete, we now exhibit results for two particular models for the initial spectrum. 



A. Conventional inflation 

We model conventional inflation as a period of power-law inflation ifllf , corresponding to an initial epoch where the 
scale factor behaves as 

a(r) = ai(rx — r)~ p ; — oo < r < t t \ . (9) 

where p > 1 (equality holds for exponential inflation) . When we write the scale factor in comoving time we get a oc t a 
with p and a related by p = a/(a — 1). 
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In this standard inflationary scenario the spectral energy density parameter for the gravitational waves can be 
calculated in the long-wavelength approximation where we assume that i;r« 1 (see the Appendix for the technical 
details) . After a straightforward but tedious calculation we obtain 
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where T is given by 



r = 



[r(- 
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sec 2 (p7r) 



When p = 1 this gives r = 1/-7T and we recover the usual result for exponential inflation. 

From the observational point of view, the extra redshift suffered by waves inside the horizon during the second stage 
of inflation will have a devastating effect on the prospects for observing them. Assuming for instance that during 
the second epoch of inflation we have exponential inflation, we will have an extra redshift of Rf , which assuming 
Ri w 10 4 gives a redshift of 10 16 . Since LIGO operates in the KHz region, which corresponds to energies of the order 
of 10 7 GeV, only when the second burst of inflation occurs at an energy scale larger than this can we ensure that 
the LIGO region will not be affected. However, in the current scenario, the second epoch of inflation will be useful 
only when it occurs at energy scales lower than this. As for LISA, depending on the scale at which the second epoch 
of inflation occurs the situation may not worsen, because LISA scales may remain superhorizon for most or all of the 
thermal inflation epoch. Nevertheless the direct detection of the gravitational wave background by LISA was already 
impossible in conventional inflationary models, because the predicted signal is too low even before thermal inflation 
does further damage. 

The exact result for the spectral energy density parameter when both the initial and second epoch of inflation 
(thermal inflation) are exponential, obtained as described in the Appendix, is plotted in Fig. [j]. This plot illustrates 
all the features discussed so far and also shows that the extra redshift introduced by thermal inflation makes the 
direct detection of the stochastic background almost impossible. In the long-wavelength approximation we get: 
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(11) 



The cosine term appearing in the third branch of the spectrum deserves some further explanation. In a full long- 
wavelength approximation we would expand the cos function keeping only the leading order term which in this case 
is 1. However, in this case it can be seen that fcr r i > 1 and the series expansion would not be appropriate. This is 
the mathematical reason behind the oscillations between the two plateaus in Fig. 

The long- wavelength approximation in Eq. (11) fits the exact result very well and it is almost indistinguishable 
from it except very close to frequencies associated with a transition and during the oscillations in the central region 
of the spectrum. The right panel in Fig. [I] shows a comparison between the exact result and the long-wavelength 
approximation in the oscillating part of the spectrum. Since the frequency of the oscillations increases with r r i <x i?ti, 
in this plot we used a small value for i? t i in order to be able to show the oscillations in detail. In the particular 
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FIG. 1. The plot on the left shows the exact spectral energy density parameter for the model of thermal inflation where 
we have exponential inflation in both inflationary eras. The spectrum for the usual model of exponential inflation without an 
extra epoch of thermal inflation is also shown for comparison purposes. The dotted lines marked LISA and LIGO II are the 
expected sensitivities of the gravitational wave detectors LISA and LIGO in the advanced configuration. Detection was already 
very problematic in conventional models of inflation Jl^], and with the inclusion of an epoch of thermal inflation the scales 
relevant for LIGO are redshifted by an additional factor Rf which makes their detection hopeless. This plot was obtained with 
pV 4 = 3 x 10 16 GeV, pf/ 4 = 10 8 GeV and i?i = I0 4 . Here, as in all the other plots in this paper, we assume h — 0.65. The 
plot on the right shows a close-up of the oscillations in the central part of the spectrum, both for the exact result and the 
long- wavelength approximation given by Eq. (|ll|). In order to be able to resolve individual oscillations, we used R tl = 30 rather 
than the realistic value of I0 4 for this illustration. All the other parameters are the same as in the left plot. We checked that 
on increasing the number of e-foldings of expansion during thermal inflation, the phase difference between the exact result and 
the approximation decreases and for R t i » 100 the two curves are almost indistinguishable. 



case shown in the figure, the oscillations in the approximation have a different phase from those in the exact result. 
However, we checked that on increasing i? t i (and therefore the frequency of the oscillations), the phase difference 
decreases, until the approximation and the exact result agree very well for values of Ra <i 100. 

The factor R^ p ^ present in the two lower frequency branches of f2 gw is due to the change in the correspondence 
between comoving and physical scales, as was discussed above. We can see that if the initial epoch is exponential 
inflation (p = 1) this effect is absent. This can be understood when we realize that the effect is a shift in the 
spectrum towards low frequencies. When the spectrum is flat, as is the case for exponential inflation, the net effect 
for frequencies which were always outside the horizon is null, while for a tilted spectrum such as the one produced by 
power-law inflation this will appear as a further redshift in the spectrum at these scales. This means that when the 
initial spectrum is generated by power-law inflation and therefore is not flat, the contribution for the anisotropies in 
the cosmic microwave background radiation will be suppressed. The contribution for the anisotropies coming from 
scalar perturbations will also be suppressed, but if the initial spectra have a different tilt there will be a net effect in 
the tensor to scalar ratio. 

One can view the effect of thermal inflation as being a 'transfer function' which processes the initial spectrum into 
a final spectrum. Its rather complicated form is given by the ratio of the final Bogolubov coefficient (3 to the initial (3 
at the beginning of the radiation epoch. In the long-wavelength approximation, the transfer function is independent 
of the initial spectrum which is being processed; for example, we expect the suppression to take exactly the same form 
in the pre big bang scenario, as we shortly illustrate. In the end the transfer function will lead to the same result as 
the full long- wavelength approximation (except possibly for a factor of 0(1)). The advantage of the transfer function 
formalism is that regardless of how the initial spectrum was produced (in this subsection conventional inflation, in 
the next the pre big bang model) for a given evolution of the Universe after the initial spectrum was produced, the 
transfer function is always the same. 

Let us now illustrate the use of the transfer function as discussed in the appendix. Neglecting constants of 0(1) the 
initial spectrum is in this case given, in terms of the number of particles Nk in mode k as defined in the appendix, as 

W '=(S^' (12) 
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where r p n = 2(r r 2 — ^ 2) +_ 7 "ri is nothing but the argument of the Hankel functions appearing in the mode functions. 
From Eqs. ( A27 ) and (|A2g|) we easily obtain 
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After all the factors involving the conformal time at the different transitions are written in terms of the parameters 
of our model with help of Eqs. (jij) to we obtain the same result as in Eq. ([Tl|), except for numerical factors. 

B. The pre big bang cosmology 

The pre big bang (PBB) is a recently-proposed cosmological scenario |||| inspired by scale factor duality one of 
the basic ideas of string theory [^3| . We will briefly summarize here some of the main features of the PBB scenario. 
For a more comprehensive review, covering both the motivation and the phenomenology, see Ref . |l4| ] . 

In the PBB scenario the Universe started its evolution from the most simple initial state conceivable, the perturbative 
vacuum of string theory. In such a state physics can be described by the tree-level low-energy effective action of string 
theory given by 



5cff 



dAd»i> ~ V((f>)] d A x (14) 



The details of the potential V(cf>) are still largely unknown, except that it goes to zero as a double exponential when 
<t> — » — 00 and it must have a non-trivial minimum at the present. 

During an initial stage, the evolution of the Universe is driven by the kinetic energy of the dilaton field. At this 
stage both the curvature scale and the dilaton increase as t — > 0_, \H\ ~ \4>\ ~ t~ x and there is in this case an exact 
solution for both the scale factor and the dilaton 

a = (-t)- 1 / 2 

(15) 

cj> = -31n(-t). 

The solution with increasing H is in fact an inflationary solution of the so-called super-inflation or pole-inflation type 
and this is one of the main peculiarities of this model. During this stage in the evolution of the Universe all the 
couplings will be very small since they are given by g = exp(0/2) sa 0. 

As soon as the curvature reaches the string scale M s = all the higher-order derivative terms in the action 

become important and the effective action Eq. (|l4j) , as well as the solution Eq. ([l5]) , are no longer valid. Unfortunately 
during the stage where the truly stringy regime is reached not much is known. What we know however is that during 
the stringy era the curvature scale must grow only mildly in order to prevent its unbounded growth with the consequent 
singularity (the Big Bang) which we are trying to avoid in this model. Therefore, whatever the exact form of the 
solution is, we expect that H w const and tj> ss const. The constraint on H will therefore lead us to some form of 
conventional inflation during the stringy regime. 

At some stage during the string phase the dilaton eventually approaches the strong coupling regime g ~ 1 and 
at this point the dilaton will freeze, leaving us with the correct value of Newton's constant, and the transition to 
conventional cosmology starting with a radiation-dominated phase will occur. 

The evolution we have described so far is based on the action Eq. (|lj) where the dilaton is non-minimally coupled 
to curvature, the so-called string frame. By applying a conformal transformation involving the dilaton we go from the 
string frame to the Einstein frame where we recover the usual Einstein-Hilbert action for gravity. When we transform 
from the string frame to the Einstein frame, the super-inflation type solution for the scale factor is transformed into 
accelerated contraction in the Einstein frame. Since physics is not changed by a conformal transformation, both 
frames are equivalent. After the string era when the dilaton freezes, the string and the Einstein frames coincide. 

One of the most remarkable predictions of the PBB model is a spectrum of gravitational waves increasing with 
frequency . The reason for this lies in the increase in the curvature during the dilaton era. Due to this characteristic 
shape, the cosmological background of gravitational waves produced in this type of model can in principle be detected 
for a broad range of parameters by the various detectors |[|. We can therefore hope that even when an epoch of 
thermal inflation is included in this model, the direct detection of the gravitational wave background may still be 
possible for certain ranges of parameters. 



G 



Without loss of generality, we will perform the entire computation of the gravitational wave background in the 
string frame, since the perturbation spectrum is frame-independent 

First of all, we must have an explicit form for the scale factor. During the string phase there is a well-known solution 
given in comoving time by Eq. (fHf), which in conformal time r takes the form a oc T ^ 1 /( l +^) and tfi oc — \/31nr. 
During the string era, however, all we know is that the curvature scale must be approximately constant. We can 
therefore model the stringy phase by an epoch of power-law inflation. The scale factor before the first radiation 
dominated era takes the form 
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where the subscript 's' denotes the beginning of the string era, a = 1/(1 + \/3), p is a free parameter and 
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For the dilaton we have 
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'0 ; T > Trl 

The constants 0du and ^ are determined by the continuity of the dilaton, while <fio is determined by the present 
value of G. We will not write them here since they will not be needed in what follows. As to (3, it is determined by 
the continuity of <f>' at the transition from the dilaton to the stringy phase and its value is given by 

V3p 3 + V3 

P = -F> = o P- ( 19 ) 

2 a 2 

This parameter will be important in the determination of the order /i of the Hankel functions appearing in the mode 
functions for the gravitons during the stringy era (see Appendix) which in this case is given by 

, = \\2p-2 P + l\ = P ^ + f 1 (20) 

After the end of the string phase, a radiation era begins and from then on the evolution of the Universe is given by 
Eq. ®. 

In this model we will have one parameter more than in the context of conventional cosmology: the total expansion 
during the string phase, R s = a rl /a s . 

In this case it is not possible to make <j>' continuous at the transition from the string phase to the radiation phase. 
This is the reason why we present the Bogolubov coefficients formalism in the Appendix in a form which does not 
require the continuity of the first derivative of R, which in this case is given by R — ae~^' 2 . 

From now on everything follows as before, except that we must use R instead of a, and the expressions for the 
Bogolubov coefficients have an extra term due to the discontinuity in cf>' . 

Before we proceed, let us just recall that in the absence of thermal inflation the spectrum for the gravitational 
waves is given by 
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Since f2 gw for modes entering the horizon after the matter-radiation equality will be far too small to have any 
observational consequences, we will forget about this part of the spectrum. The initial spectrum, i.e. the spectrum 
at the beginning of the first radiation era is 
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Applying the transfer function in Eq. ( |13| ) to this equation^ and writing all the factors of conformal time in terms of 
redshifts and energy scales we obtain: 
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where 5 = 1 + \/3 — 1 /p. Here we can see the advantage of using the transfer function approach instead of applying 
the long-wavelength approximation to the full expressions for the Bogolubov /3 coefficient: given the initial spectrum, 
for a given evolution of the universe the transfer function only needs to be calculated once. 

Figs. H and || show the exact final spectra and illustrate their dependence on the parameters of the model. In this 
model the direct detection of the stochastic background is still possible in spite of thermal inflation. As a matter of 
fact, because of the energy scale at which the second burst of inflation starts, the low-frequency peak falls exactly 
in the frequency range accessible to LISA. This is a happy coincidence and has not involved any fine tuning in the 
parameters of the model. 

The spectra show a characteristic shape with two peaks, one for frequencies around 10 2 -10 4 Hz and the other in the 
region 10 -2 -10 -4 Hz, the exact position of the peaks depending on the parameters of the model. The high-frequency 
peak was already discussed in Ref. [|l6) and it is due to the epoch of power-law inflation introduced to model the 
stringy phase. Since the part of the spectrum produced during the dilaton phase increases with frequency and the 
part of the spectrum produced in the stringy phase, which corresponds to the largest frequencies produced in this 
model, decreases with frequency, there must be a peak at the point where the two curves meet. Although the epoch 
of power-law inflation introduced to model the stringy phase is not to be taken very seriously, the fact is that during 
the string phase the curvature scale must be approximately constant, and therefore the part of the gravitational wave 
spectrum corresponding to waves produced during this period is not expected to grow with frequency. We can thus 
argue that this peak will most certainly be one of the key signatures in the gravitational wave spectrum produced by 
a PBB universe. As to the peak for frequencies around 10 -2 Hz, it is due to the era of thermal inflation. The reason 
for the oscillations is the same as in the context of conventional cosmology we discussed in the previous section. 



x We could consider the initial spectrum to be that at the beginning of the string epoch. However, in this case the transfer 
function would not be the one in Eq. ([TJ|) since it does not take into account the evolution of the gravitational waves during 
the stringy epoch. 
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FIG. 2. The spectral energy density parameter for a model of thermal inflation in the context of PBB cosmology. The plot 
on the left illustrates the dependence of f2 gw on the parameter /i, while the plot on the right shows the dependence on the 
energy scale at the end of the stringy phase. Both plots were obtained with p 1 / 4 = 10 8 GeV and Ru = 10 4 . We stress that the 
right-hand peak, near LIGO scales, violates the conditions for a perturbative calculation and is unreliable — see the text for a 
full discussion. 



From Fig. [|, we see that larger values of /i (corresponding from Eq. (|2(]) to larger values of p) correspond to larger 
values of f2 gw in the high-frequency peak. This is because the slope of logf2 gw is given by 2fi — 3. 

Increasing the energy scale at the end of the string phase shifts the entire spectrum upwards (Fig. |[ right panel). 

As to the dependence of Q gw on i? s , we see from Fig. || and the last line of Eq. (p3|) that an increase in i? s 
corresponds to a decrease in the frequency associated with the dilaton-string transition and an associated increase in 
the high-frequency branch of the spectrum which corresponds to the string era. Associated with this effect there is 
also the fact that for R s > 10 10 the amplitude of the low-frequency peak is larger than the high-frequency peak, since 
the frequency associated with the transition from thermal inflation to radiation is fixed and the branch associated 
with the dilaton era (oc / 3 ) is therefore shorter. 

Finally, the dependence of f2 gw on i? t i is somewhat unexpected. When all the other parameters are fixed, larger 
values of i?ti correspond to larger il gw in the low-frequency branch of the spectrum, contrary to what happens in 
all the other branches and what might have been expected, since the epoch of thermal inflation further redshifts the 
energy density. Once again the explanation lies in the behaviour of the frequency associated with the beginning of the 
low-frequency branch of the spectrum. Keeping all the other parameters fixed, this frequency decreases with l/i?ti 
therefore leading to a larger increase in the oscillating branch of the spectrum. 

Although the high- frequency peak also falls in the frequency range accessible to LIGO, there is a technical subtlety 
involved in the derivation of these results which may invalidate our conclusions regarding the high-frequency peak 
of the spectrum, which is that this peak falls outside the domain of application of perturbative theory as we now 
demonstrate. Our approach to the calculation of graviton production is perturbative and therefore we must ensure 
that the results obtained satisfy the domain of applicability of the perturbative calculation, which in our case means 
that we must have hij « 1. To check the validity of this condition in this particular case we need to convert our 
results for the spectral energy density parameter to the dimensionless amplitude h of the perturbations. 

We can obtain a rough estimate for h by noting that there is an equivalence between each of the two polarizations 
of the gravitational waves, + and x , and a scalar field <p: 

h +:X = VWttG ip + , x , (24) 

from which it follows that 

i fti 2 (yh) 2 \ 2 h 2 

p ^ k T^g\-o>- + —)~ mpi A^ ' (25) 

where A is a characteristic wavelength. From Eq. (B9) we can now obtain the relation we require 



fe-VAA^l (26) 
mpi / 








Given the values of f2 gw and p c at a given time we can estimate the corresponding values of the amplitude h. In the 
case under study, the problem lies in the high-frequency branch of the spectrum produced during the transition from 
the stringy era to radiation and therefore we must calculate Q gw at the beginning of the radiation era. We can obtain 
fi gw from the last line of Eq. (|23|), by taking also into account the redshift suffered by the frequency / during the 
expansion of the Universe: 

^ w d = i^(l + z cq )tt° w , (27) 

where 0^ and f2g W are the spectral energy densities of the gravitational wave background at the beginning of the 
radiation era and at the present. Using i? t j = 10 4 , if we have f2g W = 1CP 11 , the threshold for the detection of 
the high-frequency peak by LIGO, then Q, 1 ^ sw 10 10 . This is certainly a very large value and were we still in the 
regime where the perturbative calculation is valid, then this large value for would represent an inconsistency 

in our model since we are assuming f2 t ot = 1- Using pl^ 4 ~ 10 16 and noting that frequency is redshifted during the 
expansion of the Universe since the end of the stringy epoch by a factor of 10 31 , we obtain from Eq. (26) /i ra d ~ 10 9 



for frequencies of the order of kHz. This means that our results for the range of frequencies where we sec the first peak 
are not consistent with the perturbative approach and therefore they must be taken very cautiously. However, for 
the low-frequency peak generated due to the epoch of thermal inflation no such problem exists, and the perturbative 
approach employed here is justified. 

This conclusion is completely general; any stochastic background generated before thermal inflation would have to 
be non- linearly large if it were to survive thermal inflation with an amplitude detectable by LIGO. 

Notice that in this particular model, the requirement that ft < 1 at the beginning of the radiation era is much 
stronger than the constraint f2 gw < 10~ 5 imposed by nucleosynthesis |l7],[l8|]. The nucleosynthesis constraint is roughly 
equivalent to requiring h < 1 at the end of inflation for / ~ 10 10 Hz which is the high-frequency cutoff. 

As a final remark, we should notice that simultaneous detection of the stochastic gravitational wave background 
by both LISA and LIGO in the advanced phase is not ruled out. As a matter of fact the frequencies probed by both 
detectors are appropriate for the detection of both peaks in the spectrum. 



III. EARLY MATTER DOMINATION 



Another interesting possible modification to the standard cosmology is the occurrence of an early epoch of matter 
domination short after inflation. A brief early epoch of matter domination may happen during the QCD phase 
transition Q , but we will concentrate on a more prolonged epoch, such as might be caused by temporary domination 
of the Universe by moduli fields || . The scale factor after the beginning of the first radiation epoch is given by 
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a(r) = < 



a 2 T ; T ri < r < T m i 



03 (73 + t) 2 ; r m i < r < r r2 



(28) 

a.4 (T4 + r) ; r r2 < r < r oq 



^ 05 (7"5 + T 



2 



where once again the constants are chosen in order to ensure the continuity of a and a' . Notice that in this case the 
parameterization of the scale factor is slightly different from the one in the model with thermal inflation (Eq. (||)). 
Parametrizing the epoch of early matter domination by the energy scale at which it starts, H m i, and the total 
expansion of the Universe during this era i? m = a T 2/a m i, we can write all the parameters of this model in terms of 
-ff;, H m i, R m , Hq and Zgg. The value of the Hubble parameter at the time of the matter-radiation equality is in this 
case also given by Eq. (||) while for r cq we have 

r cq « {l + z cq y 1/2 H^ 1 . (29) 

Notice that although the values of r eq given by Eq. (||) for the model of thermal inflation, and by Eq. ( p9| ) for the 
model with early matter domination, are the same, the first equation is exact while the second is an approximation 
valid for z r2 ^> z cq , which is always the case in physically-realistic models. The difference between these two cases lies 
in the different parameterization for the scale factor after the beginning of the second radiation era. For the values of 
time at the end of inflation and the beginning and end of the early matter era we have 

(l + z cq ) 1/4 id /4 (l + z cq ) 1/4 id /4 ^Hl + z^B^ 

7"rl 1/0) 7ml — 1 /o 3 7V2 ^ -. / 

(HiHof 2 (H ml H ) 1/2 (H ml H Q ) 1/2 

where the approximation in r r2 is valid for R m 1. 

In the case of conventional inflation the scale factor in the initial epoch is given by Eq. (||) while in the PBB case 
the scale factor is initially given by Eq. (|l6|). 

The modifications introduced in the spectrum by the epoch of early matter domination are easy to understand in 
the light of what has been said before. First, the high-frequency cutoff will be redshifted. Since the usual matter era 
already introduces a factor (1 + Zcq) 1 / 4 in the high-frequency cutoff, we should also expect the extra redshift to be of 
the form R 



eq; 
1/4. 



R^(l + z cq )^ 



M = ,/V — ■ (31) 



As before, frequencies which are outside the horizon during the early matter domination era are not affected by it 
and we expect the spectrum to be the same for frequencies smaller than 



TJ (H m H Q ) 1/2 
hi — a r 2^r2 — —£Ti 777 (oZ) 

Zcq) 



R 3 J\l + z cq )i/± 



For frequencies in the interval / m i < / < / r i with / m i = a m i_£f m i, the spectrum will be redshifted by a factor 
R m) which corresponds to the ratio between the energy density of the gravitational waves p gw ~ a -4 and the energy 
density of matter p m ~ a -3 . The slope of this part of the spectrum is not affected since these waves have frequencies 
larger than the frequency cutoff / m i introduced by the early epoch of matter domination. For frequencies in the 
interval / r2 < / < / m i however, the early epoch of matter domination will introduce an extra factor of f~ 2 in the 
frequency dependence. This branch of the spectrum will smoothly join the redshifted part of the spectrum for high 
frequencies with the portion of the spectrum which is not affected by the early epoch of matter domination. Contrary 
to the previous case, once a given wavelength enters the horizon it will remain inside for the rest of its evolution, and 
therefore in this case there are no oscillations in the spectrum. 

Using the long-wavelength approximation we obtain 
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FIG. 4. The exact spectral energy density parameter for the model with an epoch of early matter domination in the 
context of conventional inflation (left plot) shows a suppression for / > 10 Hz, but in this case there are no oscillations. The 
suppression caused by the era of early matter domination makes the direct detection of the gravitational wave background 
almost impossible. The long-wavelength approximation is almost indistinguishable from the exact result except close to the 
transition frequencies. This plot was obtained with p\ /4 = 3 x 10 16 GcV, p l J£ = 10 10 GeV and R m = 10 4 . The right plot 
is the spectral energy density parameter for a PBB model with an epoch of early matter domination. Comparing with the 
case of thermal inflation we see the low-frequency peak does not appear in this case. Instead, the spectrum grows as / in the 
intermediate region due to the early stage of matter domination. This plot was obtained with p\( 4 = 3 x 10 16 GeV. For the 
PBB model, the same caveat that the perturbation calculation is unreliable applies as for the thermal inflation case. 
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(R m (l + z eq )) 



1/4 



(33) 



which confirms our expectations about the modifications introduced by the early epoch of matter domination. We 
also see an extra factor of R^ 1 in the two low- frequency branches of the spectrum which is, as before, due to a change 
in the correspondence between comoving and physical scales. 

Fig. ^| shows the exact result for f2 gw in a model with a short period of matter domination where the Universe started 
with an epoch of exponential inflation p — 1. We see that there is a strong suppression of modes with frequencies 
larger than / mr , the frequency corresponding to the transition from the epoch of early matter domination to the 
radiation era. The step structure is exactly what was found by Schwarz Q for the QCD transition, though much 
more prominent as our matter-dominated era is much longer. As in the previous case, the inclusion of an epoch of 
early matter domination makes the direct detection of the background of gravitational waves a very difficult, if not 
impossible, task. 

We have also analyzed a model with early matter domination in the context of PBB cosmology. Since the details 
of the model are not significantly different from those of the models previously analyzed we will omit them. The final 
result for the spectrum is shown in the Fig. ^| (right plot). In this case the low- frequency peak is absent and the early 
epoch of matter domination introduces a branch which behaves as /. Since the slope has the same sign as in the 
branch of the spectrum produced by the dilaton era, there is no peak and detection by LISA is not possible. 
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As in the previous case, the amplitude of the peak when extrapolated back in time will give rise to a value of 
/i > 1 at the beginning of the radiation phase and once again the results concerning the peak must be considered 
very cautiously. 

IV. CONCLUSIONS 

The unknown cosmological behaviour before nucleosynthesis can significantly affect the present-day spectrum of 
cosmological gravitational waves on scales probed by interferometer experiments. In this paper we have considered 
two examples of this effect, one being a short late burst of inflation and the other a prolonged period of early matter 
domination. In each case, the amplitude of gravitational waves at high frequencies is significantly suppressed compared 
to the standard cosmology. In addition, the late inflation model has a very characteristic signature, in the form of a 
series of very strong oscillations in the spectrum. 

The influence of these effects on dctectability depends on the model for the initial spectrum. In conventional 
inflation models the situation was already extremely pessimistic anyway, so arguably a further deterioration is of little 
concern. In the more speculative pre big bang models of inflation, whose ability to generate observable gravitational 
waves has been much remarked upon, the situation is more complex; while thermal inflation makes it less likely that 
gravitational waves might be seen at the LIGO frequency range, it opens up a new possibility that they may be seen at 
the LISA range. However, the situation is quite model dependent, and there is the additional worry that on LIGO scales 
the usual linear theory calculations break down, so that predictions on those scales must be treated with caution. In 
contrast, our predictions on LISA scales satisfy the linear approximation throughout. 
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APPENDIX: THE BOGOLUBOV COEFFICIENTS FORMALISM 

1. Exact results 

The mechanism responsible for the generation of the cosmological background of gravitational waves, parametric 
amplification of vacuum fluctuations, was discussed by Grishchuk in the mid seventies |ll|. To compute the spectrum 
of primordial gravitational waves we use the Bogolubov coefficient formalism [ p0| . We assume the transitions between 
the different epochs in the evolution of the Universe are instantaneous. This is in general a very good approximation 
and it will break down only for frequencies much larger than H at the instant of the transition. The accuracy of this 
approximation has been confirmed in models where the transitions are not instantaneous pl| ]. 

Since in the PBB scenario we will compute graviton production in the so-called string frame, we summarize this 
formalism in a form which is slightly more general than the one most commonly used Jl7| . The formalism laid down 
in this section can be applied to any theory whose effective action can be written in the form 

S=^ J ' J—gF^cfy^x, (Al) 

where Q is a constant, is a scalar field and TZ is the Ricci scalar. In general relativity Q = 16ttG and F = 1Z, while 
during the dilaton phase in PBB cosmology Q — l/2Ag, where A s is the string-length parameter, and F = exp(— 0)72. 
where <j> is the dilaton. We assume the geometry of the Universe is well described by a flat Friedmann-Lemaitre- 
Robertson- Walker (FLRW) line element 

ds 2 = dt 2 - a 2 (t)dx. 2 = a 2 (r) (dr 2 - dx 2 ) , (A2) 

where t is comoving time and r the conformal time obtained by dr = dt/ a. From now on we will always use conformal 
time. 

To first order, the primordial gravitational waves can be represented as a transverse-traceless perturbation h^ v of 
the FLRW metric || 
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a - o (0) 4- h 



(A3) 



where g$ is the metric of the unperturbed FLRW background. After quantization, and defining a new quantity 



R(t) = a(r) 



the hav field can be expanded as 



hii ;(x, t) = Vg'S^ f ^Yita(x, t) 



Xk,A(x,r) = a k , A 4e 4k -Vfc(r) + a kx tj < 4K ' x a4(t) , 



ik.x , 



(A4) 



(A5) 
(A6) 



where i,j — 1,2,3, the * denotes the complex conjugate, k is the comoving wave vector, k — |k| and the sum in 
Eq. (A5) is over the two polarization states described by the polarization tensor e^ . The a\\ and at (from now on 
we will drop the polarization index A) are the creation and annihilation operators for the quanta of the field h^ Vl the 
gravitons, and obey the canonical commutation relations 



Ok, »k' 



,t „t 







S(k - k') . 



(A7) 



Inserting the perturbed metric into the Einstein equations, we obtain the equation which governs the behaviour of 
the mode functions p7|,p| 



k 2 -\) in = o 



(A8) 



In most cases of interest, R oc r Q and for such an R the solution of Eq. ( A8) is given in terms of the Hankcl functions 
and H^K Imposing the requirement that in the limit r — > — oo we recover the vacuum positive frequency 



solution, /Xfc 



-ikr 



we can write [i as 



(A9) 



During its evolution, the Universe will undergo a series of transitions between different epochs with different mode 
functions, and therefore different vacuum states corresponding to each of these eras. The relation between the 
annihilation and creation operators in the initial state, and ol, and in the final state, 6k and is given by a 
Bogolubov transformation |2(J 



6k = afe ak + PI a_ 



k • 



(A10) 



where and [3k are the Bogolubov coefficients. To obtain and (3k we impose the continuity of the field hij and 
its first derivative at the instant of a transition between two different eras in the history of the Universe. In general 
relativity, where R = a, it is always possible to make R and R' continuous, but since in the context of PBB cosmology 
it may not be possible to make R' continuous, as we will see later, we assume that although R is continuous R' may be 
discontinuous at the (instantaneous) transition between two epochs. In this case, denoting with an overbar quantities 
referring to the state after the transition occurred, the Bogolubov coefficients are given by |0| 



1 



R -R 1 
R 

R' -R' ~ 
R 



(All) 
(A12) 



where all the quantities are calculated at the time of the transition. The last term in Eqs. (All) and ( A12| ) appears 
because we are not imposing the continuity of R': when R' is continuous this term is zero and we recover the usual 
expression for the Bogolubov coefficients. 

By noting that the Wronskian for the solution Eq. jASf) of Eq. ( A8 ) is given by 
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W = ii k i%-M = 2ik, (A13) 
we obtain an important relation between the Bogolubov coefficients 

\a\ 2 -\(3\ 2 = l. (A14) 
The number of particles produced in mode k is given by 

N k = \(3 k \ 2 , (A15) 
and the spectral energy density parameter is defined by 

1 

Pc dhif 



(A16) 



where / is the frequency and p c and p gw are, respectively, the critical energy density and the energy density of 
gravitational waves, and is given by 

n gw = -^^N k , (A17) 

where the angular frequency u) = 2irf = ck/a. 

When we consider two or more transitions in the history of the Universe, the final Bogolubov coefficients are given 
by the composition of the Bogolubov coefficients in each transition, according to 

a k = OL lk OL 2 k + Pik^2k > (A18) 
Pk = Pika*2 k + aikfok ■ ( Al 9) 

From now on we shall drop the index k in a and f3. Since at any instant modes with wavelengths larger than the 
horizon have not yet gone through a complete oscillation, they will not contribute to the energy density of gravitational 
waves at that given time and we are therefore forced to introduce a low-frequency cutoff for frequencies smaller than 

wif(r) = 2vriJ(r) . (A20) 

As we are assuming that all the transitions are instantaneous, a high-frequency cutoff must also be introduced. In 
an instantaneous transition modes with arbitrarily high frequency will be excited. However, a real transition will take 
a finite amount of time, and therefore these high-frequency modes are unphysical, since modes with frequencies much 
larger than the rate of expansion of the Universe will suffer many oscillations during the transition and will not be 
affected by it. The value of the high-frequency cutoff is given by p3| , p4j 

(r) = 2tt^j^. . (A21) 
a(r) 

where the subscript 'tr' refers to the time when the transition takes place. This rough argument for the introduction 
of the high-frequency cutoff can be formalized in a more rigorous way via Parker's adiabatic theorem p5|, which 



states that modes with frequencies larger than Eq. (A21) will be exponentially suppressed. Using a modification of 



the formalism presented in this section which allows for continuous transitions, and introducing a toy model for the 



transitions, it is also possible to see that modes with frequencies larger than the cutoff frequency given in Eq. (A21) 
are indeed exponentially suppressed |2lf| . 

2. Long-wavelength approximation and transfer function 

Although the application of the formalism just described is straightforward, for models involving several transitions 
the final result for f2 gw tends to become very long and not very instructive, and an approximation is due at this point. 



The high-frequency cutoff given by Eq. ( A21) can be written in comoving coordinates as 

feu = 2ira tr H u . . (A22) 

This corresponds to frequencies which are just entering the horizon at a given transition. Since k < fchf and fc^fT ~ 
(a'/a)r i=s 1, then fcr < 1 in most cases of interest and we can replace the Hankel functions appearing in the mode 
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functions (Eq. (|A9|)) by their leading terms. This approximation is usually very good for most purposes as we can 
see in Fig. [jj. Since modes satisfying fcr < 1 have wavelengths larger than the Hubble radius, this approximation is 
known as the long-wavelength approximation. 

The long-wavelength approximation can be used to generate a transfer function for the gravitational wave spectrum. 
Once we have an initial spectrum then, for a given subsequent evolution of the Universe, the same transfer function 
can always be applied to the initial spectrum, no matter what its form or origin is, to produce the spectrum at the 
present. 

Noting that the leading term in the expansion of the Hankel function is given by 



i 2 V csc(i/tt) 1 
r(i - v) (krY 



(A23) 



from Eq. ( A9) we can immediately sec that in the long wavelength approximation the mode functions are given by 

H cx (kTy u+1/2 . (A24) 

After a succession of n transitions we should expect the final (3 to be given by a product of terms of the form 
{kriY 1 ■ ■ ■ (fcT n )" n with the possibility that some of the terms may repeat them s elves in the pr oduct . After we replace 
the mode functions by their long- wavelength approximation, from Eqs. (All), ( AL2), (AlS), ( A19) we get 



P « (fcri)" 



Tn-l 



(kr n y 



"n+1 



k n+1 < k < k n 



where we dropped all the numerical factors. The index n + 1 appears in the last term since the nth transition links 
the epochs n and n + 1. The number of particles for the branch of the spectrum corresponding to the nth transition 
will therefore be given by 



N k = 



-2(1/1+1/2) -2(i/ 3 -i/ 2 ) -2(n, + i-%) u-2(vn+i+vi) u , <r k <r k 



(A25) 



From Eqs. jAll|) , (|Alg) , ( |Ajg ) we see that the final Bogolubov coefficients will be a sum of products of 

Hankel functions, each term with 2n of these, in the case of the nth transition. We can therefore ask ourselves why 
the final dependence on k is just fc-^n+i+fi) anc | d oe s not involve any of the exponents v % ■ ■ ■ v^-i which should be 
contributed by the previous transitions. The reason for this cancellation can be understood by noticing that 

H^(x) = J u {x)-iY v {x) 

where J and Y are the Bessel functions of 1st and 2nd kind. The Y function contributes the leading term x~ v 
while to lower order J ~ x v . The intermediate powers of k cancel because J ~ x v and in the end we are left with 
(3(x fc-^i+^+i). 



We can use Eq. (A25) to build a transfer function for the evolution of the gravitational waves. The initial spectrum 
(the spectrum after the first transition) takes the form 



_ T -20i+"2) £-2(i/i+i/ 2 ) 



(A26) 



We are looking for a transfer function, which we shall call G{k), such that A k = G(k)N^ 0) . From Eq. flA25| ) we can 
easily see that the required G(k) is given by: 



G(k) 



-2(, 3 -, S ),.. r -_2h-%- 1 ) T -2(, n+I -,„) r 2K +1 - 1 ,) <k< K 

-2(u 3 -vz) t _ -7 .-2(i/ n -i/„_i)^_2( 1 / n -!/ 2 ) 



n+1 


< k < 


k n 


k n 


< k < 




k 3 


< k < 




k 2 


< k < 





(A27) 



Note that G(k) is independent of the initial spectrum and can therefore, once calculated for a given cosmology, be 
applied to a variety of different choices of initial spectrum. Usually, due to requirements of continuity of the scale 
factor, the time dependence of a during the nth epoch is given in terms of r' n + r whe re r' n is a constant. As a 
consequence the argument of the Hankel functions will be k(r n + r) and therefore, in Eq. (A27) r n must be replaced 

b y T n + T n- 
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In our model of thermal inflation there is however a further complication: the cutoffs associated with the transitions 
from the first epoch of radiation to the epoch of thermal inflation, k rt i, and from thermal inflation to the final radiation 
epoch, fctir, satisfy fc rt i < fctirj contrary to the common situation where later transition will give rise to cutoffs at smaller 
frequencies. Assuming this situation involves transitions between branches i and i + 1 and between i + 1 and i + 2 
(i > 1), the transfer function will take the form 



G(k) = < 



-2(v 3 -v 2 ) -20>i-i/i_i) -2{u i+1 - Vi ) -2(v i+2 -v i+1 ) ,-2(v i+2 -v 2 ) , , , 

'2 i— 1 i i+1 "'1+2 re *~ "-1 

r -2(v 3 -v 2 ) _ _ ^-2( ! / 1 -i/i_ 1 )^-2(i/ i+ i+i/ i+ 2)^-2(i/ i +i/ 1+ i+i/ 1+2 -j/ 2 ) ki < k < k i+ i 

T -^-^... T r^-"^h- 2 ^-^ k i+1 <k< fci_! 



(A28) 



Obviously, since here we are taking only the lowest-order terms, with this approximation we cannot recover the 
oscillations seen in Figs. 0,| and |, but at least we can get the correct dependence on frequency. 



[1] K. Thorne, |gr-qc/9704042) B. Allen and J. D. Romano, Phys. Rev. D 59, 102001 (1999). 

[2] D. H. Lyth and E. D. Stewart, Phys. Rev. Lett. 75, 201 (1995); D. H. Lyth and E. D. Stewart, Phys. Rev. D 53, 784 
(1996). 

[3] A. M. Green, A. R. Liddle and A. Riotto, Phys. Rev. D 56, 7559 (1997). 

[4] D. J. Schwarz, Mod. Phys. Lett. A13, 2771 (1998). 

[5] M. Gasperini and G. Veneziano, Astropart. Phys. 1, 317 (1993). 

[6] M. Gasperini and G. Veneziano, Mod. Phys. Lett. A8, 3701 (1993). 

[7] M. Gasperini and M. Giovannini, Phys. Lett. B 282, 36 (1992); R. Brustein, M. Gasperini, M. Giovannini, and G. 
Veneziano, Phys. Lett. B 361, 45 (1995); R. Brustein, M. Gasperini, M. Giovannini, V. F. Mukhanov, and G. Veneziano, 
Phys. Rev. D 51, 6744 (1995). 
[8] M. Gasperini and M. Giovannini, Phys. Rev. D 47, 1519 (1993). 
[9] B. Allen and R. Brustein, Phys. Rev. D 55, 3260 (1997). 
[10] D. Polarski and A. A. Starobinsky, Nucl. Phys. B 385, 623 (1992). 
[11] F. Lucchin and S. Matarrese, Phys. Rev. D 32, 1316 (1985). 
[12] A. R. Liddle, Phys. Rev. D 49, 3805 (1994); (E) 51, 4603 (1995). 
[13] G. Veneziano, Phys. Lett. B 265, 287 (1991); 

[14] G. Veneziano, in Proc. of the J^th Course of the Int. School of A strophysics D. C halonge (Erice, 1995), edited by N. Sanchez 
and A. Zichichi (Kluwer Academic, Dordrecht, 1995), p. 285, hep-th/9512091; M. Gasperini, in Proc. of the 4th Course 
of the Int. School of Ast rophysics D. Cha longe (Erice, 1995), edited by N. Sanchez and A. Zichichi (Kluwer Academic, 
Dordrecht, 1995), p. 305, [hep-th/9509"l27 . 

[15] M. Gasperini and G. Veneziano, Phys. Rev. D 50, 2519 (1994). 

[16] M. Galluccio, M. Litterio and F. Occhionero, Phys. Rev. Lett. 79, 970 (1997). 

[17] M. R. de Garcia Maia, The Stochastic Background of Gravitational Waves, DPhil. Thesis (1994) University of Sussex; 

M. R. de Garcia Maia and J. D. Barrow, Phys. Rev. D 50, 6262 (1994). 
[18] B. Allen, in Relativistic Gravitation and Gravitational Radiation, edited by J.-A. Marck and J.-P. Lasota (Cambridge 

University Press, 1997). 

[19] L. P. Grishchuk, Sov. Phys. JETP 40, 409 (1975); L. P. Grishchuk, Ann. NY Acad. Sci. 302, 439 (1977). 

[20] N. D. Birrell and P. C. W. Davies, Quantum Fields in Curved Space, Cambridge Monographs on Mathematical Physics 

(Cambridge University Press, Cambridge, 1992). For applications to the case of graviton production see L. F. Abbott and 

M. B. Wise, Nucl. Phys. B244, 541 (1984) and Ref. @. 
[21] R. G. Moorhouse, A. B. Henriques, and L. E. Mendes, Phys. Rev. D 50, 2600 (1994); L. E. Mendes, A. Henriques, and 

R. G. Moorhouse, Phys. Rev. D 52, 2083 (1995), S. Koranda and B. Allen, Phys. Rev. D 52, 1902 (1995). 
[22] C. W. Misner, K. S. Thorne, and J. A. Wheeler, Gravitation (W. H. Freeman, New York, 1973). 
[23] B. Allen, Phys. Rev. D 37, 2078 (1988). 
[24] M. R. de Garcia Maia, Phys. Rev. D 48, 647 (1993). 
[25] L. Parker, Phys. Rev. 183, 1057 (1969). 



17 



